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Abstract 


A  now  method  tor  computing  the  performance  of  uniformly 
loaded  world  floors  is  presented.  The  procedure  presents  a 
floor  as  a  simple  structure  consisting  of  a  beam  supported  by 
elastic  springs.  The  method  computes  midspan  joist 
deflections  which  are  virtually  identical  to  those  obtained 
from  a  large-scale  finite  element  program,  but  at  a  fraction  of 
the  computational  effort.  Also,  computations  agree  very 
closely  with  laboratory  results.  A  simple  BASIC  program  is 
ptesented  for  implementing  the  procedure.  ^ 
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Deflections 
of  Uniformly 
Loaded  Floors 

A  Beam-Spring  Analog 

William  J.  McCutcheon,  Engineer 
Forest  Products  Laboratory,  Madison,  Wis. 


Introduction 


Wood-joist  floor  systems  have  long  been  analyzed  and  Recently,  the  FEAFLO  program  was  used  to  predict  the 

designed  by  assuming  that  the  joists  act  as  simple  beams  in  behavior  of  floors  constructed  with  joists  whose  properties 

carrying  the  design  load.  This  simple  method  neglects  many  were  determined  in  an  in-grade  survey  (Bufano  et  al.  1980; 

complex  interactions  among  a  floor’s  components— the  joists,  Vanderbilt  et  al.  1980).  These  simulations  considered  three 

sheathing,  and  connectors— which  affect  its  strength  and  stiffness  criteria;  average  midspan  joist  deflection,  soft-spot 

stiffness.  It  also  neglects  material  variability  in  that  it  deflection,  and  greatest  individual  joist  deflection  in  each 

assumes  all  joists  have  identical  mechanical  properties.  A  floor, 

simple  analytical  method  is  presented  here  which  accurately 

accounts  for  these  interactions  in  computing  floor  stiffness;  The  T-beam  model  is  computationally  much  more  efficient 

deflections  computed  by  this  method  compare  very  closely  than  a  large-scale  finite  element  program.  Because  it 

with  those  obtained  from  a  more  complex  and  costly  represents  a  floor  as  a  single  beam,  it  is  applicable  to  the 

procedure.  computation  of  average  joist  deflection;  however,  it  cannot 

directly  compute  soft-spot  or  individual  joist  deflections.  An 
Prior  research  has  addressed  the  problem  of  floor  analysis  extension  of  the  T-beam  method  to  consider  such  additional 

and  produced  analytical  methods  and  mmputer  models  which  deflection  criteria  will  provide  a  simple  and  economic 
accurately  predict  the  true  behavior  #f light-frame  floors.  alternative  for  determining  floor  performance. 

McCutcheon  (1977)  presented  an  analytical  method  for 

computing  the  stiffness  of  floors  with  partial  composite  This  paper  presents  a  method  for  predicting  the  performance 

action.  Based  on  a  T-beam  model,  the  method  accounts  for  of  wood  floors  under  uniform  loads.  Using  an  analog  which 

interlayer  slip  due  to  the  nonrigid  attachment  of  the  floor  represents  a  floor  as  a  beam  supported  by  elastic  springs,  the 

sheathing  to  the  joists.  Thompson,  Vanderbilt,  and  T-beam  method  (McCutcheon  1977)  is  extended  to  account 

Goodman  (1977)  developed  a  finite  element  computer  for  variability  of  joist  stiffness  and  two-way  action  due  to  the 

program,  FEAFLO,  for  calculating  the  performance  of  wood  cross-joist  distributional  properties  of  the  sheathing.  It  is 

floors.  Based  on  a  crossing  beam  analysis,  this  model  shown  that  this  method  computes  individual  midspan  joist 

considers  some  factors  not  included  in  the  simpler  T-beam  deflections  which  are  virtually  identical  to  those  obtained 

model,  most  notably  joist -to-joist  variability  and  two-way  from  a  large-scale  finite  element  program  (Thompson  et  al. 

action.  However,  when  the  two  methods  were  used  1977),  and  is  capable  of  accurately  predicting  the  performance 

concurrently  in  a  cooperative  study  (McCutcheon  et  al.  1981)  of  real  floors.  A  simple  BASIC  program  is  presented  for 

to  compute  the  average  deflections  of  typical  floor  implementing  the  procedure, 

configurations  subjected  to  uniform  loads,  they  produced 
results  which  were  in  very  close  agreement. 

Foschi  (1982)  presented  a  finite  element  analysis  technique 
which  includes  lateral  and  torsional  joist  deformations  as 
degrees  of  freedom  and  considers  plate  action  in  the 
sheathing;  these  are  not  considered  in  the  FEAFLO  or 
T  -beam  procedures. 


Methodology 


A  light -frame  floor  system  consists  of  multiple  parallel  wood 
joists  to  which  are  fastened  a  sheathing  made  of  plywood  or 
other  sheet  material.  In  order  to  use  the  T-beam  model  to 
account  for  composite  joist  behavior  and  to  also  consider 
two-way  action,  it  is  necessary  to  conceptually  simplify  a  floor 
in  two  steps.  At  the  first  level  of  abstraction,  the  sheathing  is 
compressed  into  a  narrow  beam  which  spans  across  the 
simply-supported  joists  (fig.  1).  This  beam  distributes  the 
load  among  the  joists.  At  the  second  level  the  joists,  which 
act  as  leaf  springs  in  supporting  the  sheathing  distribution 
beam,  are  replaced  by  simple  coil  spritigs  (fig.  2).  These 
spring  stiffnesses  account  for  joist-sheathing  composite 
action.  The  resulting  structure  is  a  beam  supported  by  elastic 
springs.  The  ends  of  the  beam  may  be  either  simply 
supported  (fig.  2a)  or  spring  supported  (fig.  2b),  depending 
upon  whether  the  end  joists  on  the  floor  are  fully  supported 
along  their  lengths  or  are  free  to  deflect,  respectively. 

Joist  “Spring”  Properties 


Floor  Width.  A 


Each  spring  in  the  analog  represents  a  floor  joist,  whose 
bending  stiffness  is  increased  by  partial  interaction  with  the 
sheathing.  The  aforementioned  T-beam  analysis 
(McCutcheon  1977)  derived  the  following  equations  to  define 
the  composite  bending  stiffness,  EE,  of  a  joist: 
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Figure  1. -Floor  system  with  sheathing 
" distribution  beam" supported  by  simply- 
supported  joists. 
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=  effective  bending  stiffness  of  the  joist, 
including  partial  composite  action 
=  bending  stiffness  of  the  joist  if  the  sheathing 
is  rigidly  attached 

=  stiffness  of  the  unconnected  joist  and 
sheathing,  taken  as  the  stiffness  of  the  bare 
joist 

=  axial  stiffness  of  the  flange  (sheathing)  and 
web  (joist)  of  the  T-beam 
=  distance  between  centroids  of  the  flange  and 
web  =  'A  (t  +  dj) 

=  sheathing  thickness 
=  joist  depth 

=  distance  between  discontinuities  (open  gaps) 
in  the  sheathing  in  the  direction  of  the  joist 
span 

=  nail  stiffness  (load/slip  ratio),  assuming 
linear  nail  behavior 
=  average  nail  spacing. 


A 


(b) 

Figure  2.  -Beam-spring  analog  structure  with  end 
joists  fully  supported  (a)  or  free  to  deflect  (b). 
(MU*  5181) 


The  midspan  deflection,  A,,  of  a  uniformly  loaded  joist  is 


Solution  of  Analog  Structure 

Analyzing  a  wood  floor  system  with  a  finite  element  program 
(5)  involves  the  solution  of  a  large  number  of  simultaneous 

equations.  A  single  joist  may  be  broken  into  as  many  as  20 
segments,  and  analysis  of  a  10-joist  floor  will  therefore 
require  the  solution  of  up  to  200  equations.  The  beam-spring 
analog  involves  just  one  degree  of  freedom  per  joist  and, 
therefore,  a  10-joist  floor  can  be  analyzed  with  just  10 
simultaneous  equations. 

The  linear  spring  constant,  k.,  required  for  the  analog  is  the  Th*  ac.tual  ec*uat'ons  wh'ch  result  from  this  method  are 

ratio  of  joist  load  to  joist  deflection:  derived  in  appendix  A,  and  a  simple  BASIC  program  for 

implementing  the  method  is  presented  in  appendix  B. 

,  El  Appendix  C  presents  complete  input  and  output  for  analyzing 

k,  =  -Zh-  =  76.8  — -  /ai  one  noor- 

'  Aj  L>  W 

Thus,  the  equivalent  spring  constant  for  each  joist  is 
computed  by  equation  (6)  after  the  joist’s  bending  stiffness, 

Eli.  is  computed  by  the  T-beam  method  summarized  by 
equations  (1-4). 

Sheathing  “Beam”  Properties 

The  bending  stiffness  of  the  analog  beam  is  equal  to  the 
stiffness  of  the  sheathing  in  the  cross-joist  direction.  If  there 
are  no  discontinuities  in  the  sheathing  in  this  direction,  that 
is,  if  the  sheathing  is  long  enough  to  extend  the  full  width  of 
the  floor,  the  analog  beam  stiffness,  El,,,  is  equal  to  X,  EjLt1. 

However,  gaps  are  usually  present  and  these  disrupt  the 
continuity  of  the  sheathing  and  reduce  its  bending  stiffness. 

The  resulting  reduced  stiffness  can  be  approximated  by 
averaging  out  the  effect  of  these  discontinuities.  For 
example,  if  the  sheathing  spans  six  joist  spaces  (as  occurs  with 
typical  16-in.  spacing  and  96-in.  lengths  of  sheathing),  a 
discontinuity  occurs  at  every  sixth  joist  in  each  row  of 
sheathing.  Thus,  on  the  average,  there  is  one-sixth  of  a  gap 
at  each  joist  crossing,  and  the  average  stiffness  is  reduced  by 
one-sixth. 

For  the  general  case,  the  bending  stiffness  may  be 
approximated  by 

EK  =  “i2~E,Lt’(1  -~F)  (7) 


A  = 

’  384  EI3 

where 

w  =  load  per  unit  length 
L  =  joist  span. 


where 

El*  =  bending  stiffness  of  analog  beam 
E,  =  bending  MOE  of  sheathing  in  cross-joist 
direction 

s  » joist  spacing 

l '  *  length  of  sheathing  (typically  96  in.) 

L,t  «  previously  defined. 


SB 


g 

IS] 


Results 


Comparison  of  Model  with 
Finite  Element  Solution 


We  used  the  finite  element  program  FEAFLO  to  determine 
the  deflection  characteristics  and  distributions  of  wood  floors 
constructed  with  joists  whose  properties  were  determined 
from  a  large-scale  sampling  program  (Bufano  et  al.  1980; 
Vanderbilt  et  al.  1980).  Each  floor  comprised  ten  2  by  8 
joists  which  were  free  to  deflect,  plus  an  additional 
continually  supported  joist  at  each  end  (fig.  2a).  Floor 
performance  was  characterized  by  three  values:  (1)  the 
average  midspan  deflection  of  the  eight  interior  joists,  (2)  the 
“soft-spot”  deflection,  defined  as  the  highest  average 
deflection  of  any  three  adjacent  joists,  and  (3)  the  maximum 
deflection  of  any  single  joist.  The  “as-graded”  joist  data 
were  used  to  analyze  138  floors  for  Douglas-fir  (green)  joists, 
and  107  floors  for  southern  pine  joists. 


From  each  group  (Douglas-fir  and  southern  pine),  1 1  floors 
were  analyzed  by  the  beam-spring  analog.  These  correspond 
to  the  stiffest  floor,  the  floors  at  the  10th  percentile  of  the 
stiffness  distribution,  the  20th  percentile, ....  90th  percentile, 
and  the  most  limber  floor. 


The  input  data  to  the  beam-spring  model  were  the  same 
values  used  in  the  FEAFLO  analyses: 


Joist  spacing,  s  =  16  inches  (in.) 

Joist  span,  L  =  157  in.  (Douglas-fir) 

L  =  154  in.  (southern  pine) 

Sheathing  thickness,  t  =  0.5782  in. 

Sheathing  axial  modulus,  along  span,  E,  =  0.9404  x  10*  lb/in. J 
Sheathing  bending  modulus,  across  span, 

Eb  =  1 .487  x  10*  lb/in.J 
Nail  stiffness  (as  assumed  in  FEAFLO  analyses), 

K,  =  30,000  Ib/in. 

Average  nail  spacing,  s„  -  6.7  in. 

Joist  size,  nominal  2  by  8,  actual  dimensions  input  for  each 
joist 

Joist  modulus  input  for  each  joist 
Length  of  sheathing,  l'  =  9 6  in. 


The  T-beam  model  can  account  for  the  presence  of  complete 
discontinuities,  i.e.,  open  gaps,  in  the  sheathing;  but  the 
FEAFLO  analyses  assumed  flexible  gaps  between  the  48-inch- 
wide  pieces  of  plywood.  Therefore,  it  was  necessary  to 
determine  empirically  an  equivalent  distance  (greater  than  48 
in.)  between  open  gaps,  L Previously  McCutcheon  et  al. 
(1981)  determined  that  L '  =  96  inches  can  be  used  to 
account  for  flexible  gaps  with  a  stiffness  of  10,000  lb/in.  Vin. 
In  these  analyses,  a  gap  stiffness  of  5,000  lb/in.  Vin.  was  used 
in  the  FEAFLO  input,  and  L '  =  75  inches  gave  good 
agreement  between  the  two  methods. 


The  FEAFLO  and  beam-spring  analyses  (table  1)  give 
virtually  identical  results  for  all  three  deflection  criteria- 
average  interior  deflection,  soft  spot,  and  maximum  joist. 
Complete  FEAFLO  data  were  available  for  southern  pine 
floor  No.  74,  and  are  compared  to  the  beam-spring  data  in 
table  2.  These  results  are  presented  graphically  in  figure  3, 
which  also  shows  the  results  of  a  traditional  bare  joist 
analysis,  which  assumes  no  composite  action  and  no  two-way 
action. 
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Figure  3. -Comparison  of  bare  joist,  finite  element 
(FEAFLO),  and  beam-spring  model  analyses  of 
southern  pine  floor  No.  74.  (ML84  5182) 
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Tabic  1 . -Comparison  of  floor  deflections  as  computed  by  finite  element  program  (FEAFLO)  and  beam-spring  analog 


Floor  No. 

Average  eight  interior  joists 

Soft  spot  (three  joist  average) 

Maximum  (at  joist  number) 

FEAFLO 

Beam-spring 

Difference 

FEAFLO 

Beam-spring 

Difference 

FEAFLO 

Beam-spring 

Difference 

In. 

In. 

Pet 

In. 

In. 

Pci 

In. 

In. 

Pci 

DOUGLAS-F1R 

82 

+  0.7 

0.2673 

0.2743 

+  2.6 

0.3078  (2) 

0.3142  (2) 

+  2.1 

69 

.2989 

+  0.6 

.3119 

+  0.3 

.3220  (8) 

.3242  (8) 

+  0.7 

.3111 

+  0.6 

.3189 

.3211 

+  0.7 

.3273  (7) 

.3285  (7) 

+  0.4 

15 

.3189 

.3209 

+  0.6 

.3355 

.3383 

+  0.8 

.3441  (3) 

.3450  (3) 

+  0.3 

111 

.3296 

+  0.4 

.3631 

.3627 

-0.1 

.3755  (8) 

.3787  (8) 

+  0.9 

99 

.3421 

.3443 

+  0.6 

.3521 

.3554 

+  0.9 

.3580  (3) 

.3601  (3) 

+  0.6 

.3483 

+  0.6 

.3712 

.3721 

+  0.2 

.3758  (5) 

.3749  (8) 

-0.2 

29 

.3576 

.3586 

+  0.3 

.3730 

.3720 

-0.3 

.3996  (9) 

.3862  (9) 

-3.4 

55 

.3735 

.3756 

+  0.6 

.4227 

+  0.4 

.4460(3) 

.4401  (3) 

-1.3 

62 

.4047 

+  0.3 

.4332 

.4355 

+  0.5 

.4410  (8) 

.4489  (8) 

+  1.8 

35 

.4591 

+  0.4 

.4979 

.4990 

+  0.2 

.5174  (4) 

.5146  (4) 

SOUTHERN  PINE 

35 

.2228 

+  0.5 

.2366 

.2390 

+  1.0 

.2545  (2) 

.2593  (2) 

+  1.9 

91 

.2814 

.2797 

-0.6 

.2977 

.3022 

+  1.5 

.3419  (9) 

.3218  (8) 

-5.9 

37 

.2914 

+  0.4 

.3102 

-0.1 

.3220  (7) 

.3224  (7) 

jMtfirlBil 

+  1.4 

.3166 

.3140 

-0.8 

.3438  (9) 

.3301  (3) 

36 

.3169 

+  0.3 

.3472 

.3484 

+  0.3 

.3628  (6) 

.3572  (6) 

-1.5 

24 

.3244 

.3247 

+  0.1 

.3675 

.3658 

-0.5 

.3822  (9) 

.3810 (8) 

'74 

.3377 

.3383 

-0.2 

.3713 

.3708 

-0.1 

.3964  (9) 

.3842  (9) 

-3.1 

29 

.3465 

.3474 

+  0.3 

.3653 

.3643 

-0.3 

.3832  (9) 

.3770  (8) 

-1.6 

31 

.3545 

.3549 

+  0.1 

.3939 

.3912 

-0.5 

.3993  (9) 

.3989  (8) 

-0.1 

12 

.3792 

.3792 

.4120 

.4097 

-0.6 

.4700  (9) 

.4508  (9) 

-4.1 

65 

.4494 

.4512 

+  0.4 

.4852 

.4876 

+  0.5 

.4993  (3) 

.4994  (3) 

'Complete  analyses  of  this  floor  are  presented  in  table  2  and  Figure  3. 


Table  2. -Complete  analysis  of  southern  pine  floor  No.  74 


Joist 

No. 


_ Joist  properties _ 

Width  Depth  Modulus  of 

elasticity 


Joist  deflection 

Beam-spring  Difference 


Bare  joist 


FEAFLO 


Comparison  of  Model  with  Test  Data 

In  developing  the  original  T-beam  analysis  method 
(McCutcheon  1977),  the  performance  of  seven  floors  was 
evaluated  experimentally.  Of  these  seven,  two  floors  were 
intentionally  constructed  with  joists  which  had  high  degrees  of 
variability  in  their  stiffnesses.  (The  other  five  had  nearly 
uniform  joist  properties.)  Designated  N-3  for  the  floor  with 
nailed  sheathing  and  G-3  for  the  floor  with  the  sheathing 
attached  by  means  of  a  rigid  adhesive,  these  two  are  of  the 
greatest  interest  because  they  are  most  useful  in  assessing  the 
beam-spring  model's  ability  to  properly  account  for  joist -to- 
joist  variability. 

The  properties  of  the  floors  and  the  corresponding  input  data 
to  the  model  were: 

Number  of  joists,  all  free  to  deflect  (fig.  2b),  9 

Joist  span,  L  =  144  in. 

Joist  spacing,  s  =  15.8125  in.  (average  of  six  spaces  at  16 
in.  and  two  at  15.25  in.) 

Joist  size,  2  by  8  in.  nominal,  1.5  by  7.25  in.  actual 

Joist  moduli,  determined  individually 

Sheathing  thickness,  t  =  0.625  in. 

Sheathing  properties,  determined  experimentally: 
axial  modulus,  along  span,  E,  =  0.5  x  10*  lb/in. 1 
bending  modulus,  across  span,  Eb  =  1  x  104  lb/in. 1 

Length  of  sheathing,  l '  =  %  in. 

Width  of  sheathing.  L  *  =  48  in.  (N-3) 

=  144  in.  (G-3,  glued  tongue-and- 
groove  edges) 

(N-3)  nail  stiffness  (as  computed  by  McCutcheon  1977), 
k„  =  9,400  lb/in.  average  nail  spacing,  s„  =  7.43  in. 

(G-3)  rigid  adhesive,  kn  =  1  x  10”  lb/in. 

s„  =  1 

Except  for  end  joists  I  and  9,  the  analog  was  an  excellent 
predictor  of  floor  N-3’s  performance,  especially  at  the  higher 
load  (table  3,  fig.  4).  However,  it  is  unlikely  that  the  same 
nail-slip  modulus  would  apply  to  both  load  levels.  Use  of  a 
higher  modulus  at  the  lower  load,  to  account  for  the  steeper 
slope  of  a  curvilinear  load-slip  curve,  would  bring  the 
computed  results  into  even  closer  agreement  with  the 
experimental. 

The  analog  also  did  a  good  job  of  predicting  the  glued  floor’s 
performance  (table  4,  fig.  5),  except  for  the  anomaly  at  joists 
7  and  8,  where  the  experimental  deflection  of  the  apparently 
stiffer  joist,  No.  8,  was  greater  than  that  of  the  more  limber 
joist.  No.  7. 


JOIST  NO. 


Figure  4. -Experimental  and  analytical  midspan 
joist  deflections  for  floor  N-3.  (ML84  5183 ) 


JOIST  NO. 


Figure  5.  -Experimental  and  analytical  midspan 
joist  deflections  for  floor  G-3.  CDeflection  gage 
No.  3  malfunctioned)  (ML84  5184) 
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-Comparison  of  experimental  and  analytical  performance  of  floor  N-3 


Joist  deflection  at 
49.2  pounds  per 
square  foot 


Joist  deflection  at 
98.8  pounds  per 
square  foot 


Differ¬ 
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Beam¬ 
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Pel  , 

-24.0 

0.386 

0.279 

-27.7  « 

+  5.6 

.639 

.644 

+  0.8  ' 

+  6.0 

.734 

.742 

+  1.1  \ 

+  7.1 

.646 

.639 

-1.1  ; 

+  9.3 

.638 

.636 

-0.3  ] 

+  2.4 

.692 

.680 

-1.7 

+  2.4 

.594 

.599 

-0.8 

+  3.4 

.483 

.483 

0 

.381 

.304 

-20.2 

Table  4. -Comparison  of  experimental  and  analytical  performance  of  floor  G-3 


Discussion 


Conclusion 


The  input  to  this  method  consists  of  the  physical  dimensions 
and  mechanical  properties  of  the  floor's  components. 

Joist  bending  stiffnesses  can  be  easily  determined  from  simple 
beam  flexure  tests,  such  as  those  specified  by  ASTM  Standard 
D  198  (ASTM  1976).  Similarly,  sheathing  properties  (axial 
stiffness  parallel  to  the  joists  and  bending  stiffness 
perpendicular)  can  be  measured  from  axial  and  bending  tests 
of  appropriately  oriented  sheathing  specimens.  Standards  for 
these  tests  may  be  found  in  ASTM  D  3500,  D  3501,  and 
D  3043. 

Appropriate  nail  stiffnesses  may  be  more  difficult  to  define. 
The  performance  of  mechanical  fasteners  is  the  subject  of 
continuing  research  by  a  number  of  investigators.  As 
research  in  this  area  progresses,  it  should  be  possible  to 
determine  which  test  procedures  are  most  appropriate  for 
defining  nail  performance  in  various  structural  applications. 
The  current  state  of  the  art  for  testing  mechanical  fasteners  is 
given  in  ASTM  D  1761. 

However,  a  floor’s  performance  is  relatively  insensitive  to 
changes  in  the  interlayer  fastener  stiffness.  For  example,  if 
the  nail  stiffness  k„  in  experimental  floor  N-3  (table  3)  is 
increased  by  a  factor  of  2.5  from  9,400  to  23,500,  the 
computed  joist  deflections  will  decrease  by  about  6  percent. 
Similarly,  if  k„  in  southern  pine  floor  No.  74  (table  2)  is 
reduced  by  2.5  from  30,000  to  12,000,  the  deflections  will 
increase  by  an  average  of  only  16  percent. 


A  simple  method  has  been  presented  for  computing  the 
midspan  joist  deflections  in  a  uniformly  loaded  wood  joist 
floor.  A  T-beam  model  accounts  for  the  composite  action 
between  the  joists  and  sheathing,  and  a  beam-on-elastic - 
springs  model  accounts  for  joist  variability  and  two-way 
action. 

The  method  yields  calculated  results  which  are  virtually 
identical  to  those  obtained  from  a  large-scale  finite  element 
analysis,  but  at  a  fraction  of  the  computational  effort,  and 
results  also  compared  closely  with  experimental  data. 
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Appendix  A 

Derivation  of  Beam-Spring  Equations 


Simply  Supported  Beam 


A  simply  supported  beam  on  elastic  springs  is  shown  in  figure 
A-l.  It'  A*  is  the  deflection  at  node  i  due  to  the  uniform  load 
w,  and  A,.  is  the  deflection  due  to  spring  force  j,  then  the 
deflection'at  spring  i.  A,,  is: 


A  =  A^  +  I  A,  +  A,.  +  2  A,. 

j=l  '  '  j=i+ 1  1 


-  I  \  +  A  -  A  -  2  \  =  A, 

j=l  1  '  j=i+l  1 


Expanding  the  expressions: 


i  -  1  Fa  (/  -  a,) 

2 - (2a  /  -  a  2  -  a,1)  +  A, 

.  .  1  1  )  1 


F,a,2(/  -  a,)’  n  F  (l  -  a  )a, 

- -  i  -  (2a  f  -  a  2  -  a  2) 

3E1/  j=i+1  6EI/  1  1  ’ 


-(!'  -  2/a,1  +  a-’) 


a,  =  is 

a,  =  js 

I  =  (n  +  1  )s 

F,  =  -  k,A, 


and,  defining 


El  El 

k'  =  =  (77 


(n  +  1) V 


Equation  (A-3)  becomes: 


I  /3A,j(l  -  i)(2i  -  i  -  j)  +  2A,  [3  +  ft,  id  -  i)] 
i=  i 


♦  *  P,  A/l  —  j )  i  ( 2  j  -  i  -  j  ) 

i  =  i  +  1 
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K,>  Kz>  K,; 


t  I  I 


Fl  F2  F3 


(n  +  l)  ®  S=  & 


Figure  A-l. -Beam  on  elastic  springs  with  simply- 
supported  ends.  (ML84  5178) 


Thus,  from  equation  (A-4),  the  deflections  of  the  beam 
(floor)  and  every  spring  (joist)  can  be  determined  by  solving  n 
simultaneous  equations  of  the  form: 


[K1  {A}  =  {W} 


A  A  AAA 

K„  =  ftj  (1  -  0(2i  -  i  -  j  ) 


for  i  >  j 


=  2[3  +  /3j  i  (1  -i)] 
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Beam  with  “Free”  Ends 

A  beam  with  its  ends  spring-supported  is  shown  in  figure  A-2. 
The  derivation  for  this  case  is  similar  to  that  for  the  simply- 
supported  beam  except  that  additional  terms  are  required  to 
account  for  the  movement  of  the  beam  ends.  The  deflection 
at  node  i  is 


i-1  n- 1 

ft-A»+IV+V  +  1  V  + 

j  =2  1  1  j  =  i  + 1  1 


w  ~  »■> 
2  *■  k, 
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(A-5) 


Expanding  the  expressions  for  AK  and  substituting: 

a,  =  (i  -  l)s 
a,  =  0  -  Ds 
I  =  (n  -  l)s 
Fj  =-kA 
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Again,  the  problem  reduces  to  a  set  of  n  simultaneous 
equations  of  the  form 

[K]  {A}  =  {W} 
where 

Kit  =  K,„  =  0 

except 

Km  =  K„„  =  6 


Figure  A-2. -Beam  on  elastic  springs  with  "free" 
ends.  (ML84  5179) 
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Appendix  B: 

Basic  Program  for  Floor  Analysis 


The  following  BASIC  program  implements  the  method 
presented.  It  is  written  to  run  on  a  Sharp  model  1500  pocket 
computer  but  can  easily  be  modified  for  other  machines.  The 
program  prompts  for  input  data  and  requires  that  all  units  be 
consistent;  e.g.,  if  span  is  in  inches,  then  joist  dimensions 
must  also  be  in  inches,  and  if  joist  E  is  in  pounds  per  square 
inch,  then  uniform  load  must  also  be  in  pounds  per  square 
inch,  not  pounds  per  square  foot. 

The  program  is  organized  as  follows: 

Segment  “F”,  statements  10-40 
Initiate  new  problem  (clear  memory). 

Input  floor  type  (simply  supported  or  “free”). 

Input  number  of  joists. 

“J”.  100-170 

Input  joist  data  (spacing,  span,  widths,  depths,  elastic 
moduli). 

Note:  Input  zero  for  properties  which  are  not  identical  for 
all  joists. 

“S”,  200-250 

Input  sheathing  data  (thickness,  Eaxial,  E^,^,  gap 
spacings). 


“N*\  300-320 

Input  nail  data  (stiffness,  spacing). 

“D”,  400-790 
Solution: 

400-450  computes  T-beam  stiffnesses. 

455  computes  load  vector. 

460-495  assembles  stiffness  matrix. 

500-507  modifies  stiffness  matrix  for  “free”  ends. 

5 10-740  solves  simultaneous  equations. 

750-790  computes  average,  soft  spot,  and  maximum 
deflection  coefficients. 

“L”.  795 

Input  uniform  load. 

“A”,  800-860 
Display  joist  deflections. 

999  End. 


Program  Listing 


10:  “F”:CLEAR  :WA1T  0 
15:  CLS  :PR1NT  “FLOOR  TYPE(SS  or  FRee) 
”;:INPUT  F$ 

16:  IF  F$  =  “SS”GOTO  30 
17:  IF  F$  =  “FR”GOTO  30 
20:  GOTO  15 

30:  CLS  : PRINT  “No.  of  JOISTS  =”;:INPUT  N 
40:  DIM  B(N),D(N),E(N),K(N,N),Q(N),BT(N) 

100:  “J”: J$  =  “JOIST”:WAIT  0 

102:  CLS  :PRINT  J$  +  “SPACING  =  **;:INPUT  S 

103:  LL  =  (N-  1  +(F$  =  “SS”)*2)*S 

105:  CLS  : PRINT  J$+  “  SPAN  =  ”;:INPUT  L 

110:  CLS  :PRINT  JS  +  “  WIDTH  (all)  =  *’;:INPUT  B 

115:  CLS  : PRINT  J$  +  “  DEPTH  (all)  =  ”;:INPUT  D 

120:  CLS  : PRINT  JS  +  “  MoE  (all)  =  ”;:INPUT  E 

125:  USING”###” 

130:  FOR  J  =  1TO  N 

140:  IF  B>0LET  B(J)  =B:GOTO  150 

145:  CLS  : PRINT  J$;J;“  WIDTH  =  ”;:INPUT  B(J) 

150:  IF  D>0LET  D(J)  =  D:GOTO  160 

155:  CLS  :PRINT  J$;J;“  DEPTH  =  ”;:INPUT  D(J) 

160:  IF  E>0LET  E(J)  =  E:GOTO  170 

165:  CLS  :PRINT  J$;J;“  MoE  =  ”;:INPUT  E(J) 

170:  NEXT  J 

200:  “S”:S$  =  “SHTHNG”:WAIT  0 

210:  CLS  :PRINT  S$  +  “THICKNESS  =  ”, -.INPUT  T 

220:  CLS  :PRINT  SS+  “E(axial)  =  ”;:INPUT  EA 

230:  CLS  :PRINT  S$  +  “E(bndg)  =  ’’"INPUT  EB 

240:  CLS  :PRINT  “GAP(along  joist)  =  ”;:INPUT  LJ 

250:  CLS  :PRINT  “GAP(perp  to  joist)  =  ”;:INPUT  LS 

300:  “N”:NS  =  “NAIL  ”:WAIT  0 

310:  CLS  :PRINT  N$+  “STIFFNESS  =  ”;:INPUT  KN 

320:  CLS  :PRINT  N$  +  “SPACING  =  ”;:INPUT  SN 

400:  “D”:IS  =  EB*T  A3/12*L*(I  -S/LS) 

405:  FOR  J  =  1TO  N: JJ  =  (J  -  1  +  (FS  =  “SS”))trS/LL 
410:  A1  =  EA  *T  *S:  A2  =  E(  J)  *B( J)  *D(  J):  IU  =  A2 * 
D(J)A2/12 

415:  IF  FS  =  “FR”AND(J  =  1  OR  J  =  N)LET  A1  =  Al/2 
420:  H2  =  (D(J)  +  T)  A  2/4:  IR  =  IU  +  A 1  *A2/(A  1  +  A2)  *H2 
430:  FD  =  10/(H2*KN/SN*LJ  A  2*IR/IU/(IR  -  IU)  +  10) 
440:  EI  =  IR/(1+FD*(IR/IU-1)) 

450:  BT  =  76.8*EI/IS*(LL/L)  A  3:BT(J)  =  BT 
455:  Q(J)  =  LL  A  4*L/4/IS*JJ*(l  -  2*JJ  A  2  +  JJ  A  3) 

460:  FOR  I  =  1TO  N:II  =  (I  -  1  +  (FS  =  “SS”))*S/LL 
470:  IF  I  >  JLET  K(I,J)  =  BT*JJ*(1  -II)* 

(2*11  -  II A  2- JJ  A  2) 

480:  IF  I  =  JLET  K(I,J)  =  2*(3  +  BT*II  A  2*(1  —  II)  A  2) 
490:  IF  I <  JLET  K(I,J)  =  BT*(1  -  JJ)*II* 

(2* JJ  —  II A  2- JJa  2) 

495:  NEXT  I:NEXT  J 
500:  IF  FS  =  “SS”GOTO  510 


501:  FOR  1  =  1TO  N:II  =  (I  —  1)/(N-  1) 

502:  Q(I)  =  Q(I)  +  3*LL  A  4*L/IS*((1  -  II)/BT(1)  + 
II/BT(N)) 

503:  FOR  J  =  1TO  N:JJ  =  (J  -  1  )/(N  -  1) 

504:  K(I,J)  =  K(I,J)  +  6*BT(J)*((1  -  II)*(1  -JJ)/BT 
(1)  +  II*JJ/BT(N)) 

505:  IF  J  =  lOR  J  =  N  LET  K(I, J)  =  0 

506:  NEXT  J:NEXT  I 

507:  K(1,1)  =  6:K(N,N)  =  6 

510:  FOR  K=  1TO  N-  1 

520:  FOR  I  =  K  +  1TO  N 

530:  IF  K  =  1GOTO  570 

540:  FOR  J  =  1TO  K-  1 

550:  K(1,K)  =  K(I,K)  -  K(I,J)*K(J,K) 

560:  NEXT  J 

570:  K(I,K)  =  K(I,K)/K(K,K) 

580:  NEXT  I 

590:  FOR  I  =  K+lTO  N 

600:  FOR  J  =  1TO  K 

610:  K(K  +  1,1)  =  K(K  +  1,1)  -  K(K  +  1,J)*K(J,I) 

620:  NEXT  J:NEXT  I 
630:  NEXT  K 
640:  FOR  I  =  2TO  N 
650:  FOR  K  =  1  TO  1-1 
660:  Q(I)  =  Q(I)-K(I,K)*Q(K) 

670:  NEXT  K:NEXT  I 
680:  FOR  I  =  NTO  1  STEP  -1 
690:  IF  I  =  NGOTO  730 
700:  FOR  K  =  I  +  1TO  N 
710:  Q(I)  =  Q(I)  -  K(I,K)*Q(K) 

720:  NEXT  K 
730:  Q(I)  =  Q(I)/K(I,I) 

740:  NEXT  I 
750:  MX  =  0:SS  =  0:SI  =  0 
760:  FOR  J  =  2TO  N  -  1 
770:  IF  Q( J)  >  MXLET  MX  =  Q( J): AM  =  J 
780:  X  =  (Q(J  - 1)  +  Q(J)  +  Q(J  +  1))/3:IFX>SSLET 
SS  =  X:AS  =  J 

790:  SI  =  SI  +  Q(J):NEXT  J:BEEP  I 
795:  "L”:WAIT  0:CLS:PRINT  “UNIFORM  LOAD  = 
”;:INPUT  Q 

800:  "A”:US  =  “#  #.#  #  ft  #”:WAIT  :CLS 

810:  PRINT  “AVE.  DEFL  =”;USING  U$;SI/(N - 2)*Q 

820:  PRINT  “SOFT  SPOT  =  ”;USING  U$;SS*Q; 

“  @  ”;USING  ;AS 

830:  PRINT  “MAX.  DEFL  =  ”;USING  U$;MX*Q; 

"  @  ”;USING  ;AM 
840:  FOR  I  =  1TO  N 

850:  PRINT  “DEFL”;  USING  “#  #  #”;1;  “  =”;US1NG 
US;Q(I)*Q 
860:  NEXT  I 
999:  END 
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Appendix  C: 

Sample  Input  and  Output 


To  demonstrate  the  use  of  the  BASIC  program,  complete 
input  and  output  are  presented  for  experimental  floor  N-3 
(table  3). 


Program  Prompt 

Input 

Output 

FLOOR  TYPE  (SS  OR  FR)? 

FR 

AVE.  DEFL  =  0.3145 

NO.  OF  JOISTS  =  ? 

9 

SOFT  SPOT  =  0.3360  @  3 

JOIST  SPACING  =  ? 

15.8125 

MAX.  DEFL  =  0.3696  @  3 

JOIST  SPAN  =  ? 

144 

DEFL  1  =  0.1388 

JOIST  WIDTH  (all)  =  ?' 

1.5 

DEFL  2  =  0.3204 

JOIST  DEPTH  (all)  =  ?' 

7.25 

DEFL  3  =  0.3696 

JOIST  MOE  (all)  =  ?* 

0 

DEFL  4  =  0.3181 

JOIST  1  MOE  =  ? 

2.32E6 

DEFL  5  =  0.3167 

JOIST  2  MOE  =  ? 

1.85E6 

DEFL  6  =  0.3385 

JOIST  3  MOE  =  ? 

1.29E6 

DEFL  7  =  0.2981 

JOIST  4  MOE  =  ? 

2.40E6 

DEFL  8  =  0.2403 

JOIST  5  MOE  =  ? 

2.03E6 

DEFL  9  =  0.1511 

JOIST  6  MOE  =  ? 

1.42E6 

JOIST  7  MOE  =  ? 

2.27E6 

JOIST  8  MOE  =  ? 

2.67E6 

JOIST  9  MOE  =  ? 

1.74E6 

SHTHNG  THICKNESS  =  ? 

.625 

SHTHNG  E  (axial)  =  ? 

.5E6 

SHTHNG  E  (bndg)  =  ? 

1E6 

GAP  (along  joist)  =  ? 

48 

GAP  (perp  to  joist)  =  ? 

96 

NAIL  STIFFNESS  =  ? 

9400 

NAIL  SPACING  =  ? 

7.43 

(Program  computes  for  approximately  2  minutes) 
UNIFORM  LOAD  =  ?  49.2/144 


'Enter  zero  if  joist  properties  are  not  all  identical;  in  this  case, 
program  prompts  for  individual  entry. 


kV 


V., 


ua 


Appendix  D: 

'-'V7v' 

i 

£ 

Program  Variables 

Prof  nun 

Text 

Si 

Variable 

Variable 

Meaning 

y 

FS 

Floor  type:  SS  for  simply  supported 

FR  for  “free” 

* 

N 

n 

Number  of  joists 

B,B(J) 

Joist  widths 

D,D(J) 

ft 

Joist  depths 

v: 

E,E(J) 

Joist  modulus  of  elasticity 

K(I,J) 

Kii 

Stiffness  matrix  element 

Q(J) 

w, 

Load  vector  element  (solution) 

;/•; 

Q(J) 

ft 

Displacements  (output) 

\'v 

Q 

W 

Uniform  load 

» ;j 

BT,BT(J) 

ft 

Nondimensional  spring  stiffness 

i 

S 

s 

Joist  spacing 

Xr 

LL 

1 

Floor  width 

K 

L 

L 

Floor  span 

>} 

T 

t 

Sheathing  thickness 

EA 

E. 

Sheathing  axial  modulus  of  elasticity, 
along  joist 

■v 

jj 

EB 

Efc 

Sheathing  bending  modulus  of  elasticity, 
across  joist 

LJ 

L' 

Sheathing  gap,  along  joist 

LS 

I ' 

Sheathing  gap,  across  joist 

KN 

K 

Nail  stiffness 

'•'5 

SN 

S„ 

Nail  spacing 

IS 

EIb,EI 

Sheathing  El 

■> 

II 

1 

Nondimensional  joist  location 

pu 

JJ 

J 

Nondimensional  joist  location 

A1 

(EA), 

Sheathing  axial  stiffness 

K 

c* 

A2 

(EA)2 

Joist  axial  stiffness 

H2 

h1 

Square  of  distance  between  centroids  of 

T-beam  flange  and  web 

►.wj 

& 

IU 

Ely 

Stiffness  of  bare  joist 

IR 

eir 

Stiffness  of  rigidly  connected  T-beam 

El 

EIi 

Joist  T-beam  stiffness 

*i 

FD 

See  equation  (3) 

MX 

Maximum  joist  deflection  coefficient 

AM 

Location  ( joist  No.)  of  maximum  deflection 

SS 

“Soft-spot”  deflection  coefficient 

AS 

Location  (central  joist  No.)  of  soft  spot 

\v 

SI 

Sum  of  interior  joist  deflection  coefficients 

X 

Average  deflection  coefficient,  three  joists 
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